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We explore the influence of the deformation on the nuclear matrix elements of the neutrinoless 
double beta decay (NME), concluding that the difference in deformation -or more generally in the 
amount of quadrupole correlations- between parent and grand daughter nuclei quenches strongly the 
decay. We correlate these differences with the seniority structure of the nuclear wave functions. In 
this context, we examine the present discrepancies between the NME's obtained in the framework 
of the Interacting Shell Model and the Quasiparticle RPA. In our view, part of the discrepancy can 
be due to the limitations of the spherical QRPA in treating nuclei which have strong quadrupole 
correlations. We surmise that the NME's in a basis of generalized seniority are approximately model 
independent, i. e. they are "universal". 
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I. INTRODUCTION 



The double beta decay is a rare weak process which takes place between two even-even isobars when the single 
beta decay is energetically forbidden or hindered by large spin difference. The two neutrino beta decay is a second 
order weak process — the reason of its low rate — , and has been measured in a few nuclei. The 0vf]/3 decay is analog 
but requires neutrinos to be Majorana fermions. With the exception of one unconfirmed claim it has never been 
observed, and currently there is a number of experiments either taking place or expected for the near future — see 
e.g. ref. — devoted to detect this process and to set up firmly the nature of neutrinos. Furthermore, the OvjSp 
decay is also sensitive to the absolute scale of the neutrino mass, and hence to the mass hierarchy. Since the half-life 
of the decay is determined, together with the masses, by the nuclear matrix element for the process, its knowledge is 
essential to predict the most favorable decays and, once detection is achieved, to settle the neutrino mass scale and 
hierarchy. 

Two different methods were traditionally used to calculate the NME's for Qvfip decays, the quasiparticle random- 
phase approximation and the shelljnodel in large valence spaces (ISM). The QRPA has produced results for most of 
the possible emitters since long d-tsj. The ISM, that was limited to a few cases till recently Q, can nowadays describe 
(or will do it shortly) all the experimentally relevant decays but one, the decay of ^^°Nd. Other approaches, that 
share a common prescription for the transition operator (including higher order corrections), and for the treatment of 
the short range correlations (SRC) and the finite size effects, are the Interacting Boson Model Q, and the Projected 
Hartree Fock Bogolyuvov method ^] . 

The expression for the half- life of the Oi^/3/3 decay can be written as 

(T°;f (0+ ^ 0+))"' = Goi \M^^^^\' (^)' ' (1) 

where (w^^) = C^g^.TO/c| is the effective neutrino mass, a combination of the neutrino mass eigenvalues mfe. U 
is the neutrino mixing matrix and Gqi is a kinematic factor dependent on the charge, mass and available energy of 
the process. M^'^^^ is the nuclear matrix element of the neutrinoless double beta decay operator, which has Fermi, 
Gamow- Teller and Tensor components. The kinematic factor Gqi depends on the value of the coupling constant qa- 
In addition, some calculations use different values of ro in the formula R=ro A^/'^. It is therefore convenient to define: 



M' '^'^PP = ' ^ M0-/3^ (2) 

In this way the theoretical M'^'^'^^'s are directly comparable among them irrespective of the values of 17^1 and 
ro employed in their calculation, since they share a common Goi factor — the one computed with = 1.25 and 
ro=1.2 fm. Thus, the translation of the M' o^^/^'s into half-lives is transparent. 



II. PAIRING AND QUADRUPOLE; THE INFLUENCE OF DEFORMATION 



An important issue regarding the Oi^/3/3 decay is the role of the correlations; pairing that drives the nucleus toward 
a superfluid state and quadrupole that favors deformed intrinsic shapes. It has been show recently that the 2^/3/3 
is hindered by the difference in deformation between the initial and final nuclei [l3j[il|- For the neutrinoless mode, 
the calculations Q indicate that the pairing interaction favors the decay and that, consequently, the truncations in 
seniority, which quench the pair breaking action of the quadrupole correlations, produce an overestimation of the 
values of the NME's. On the other hand, the NME's are also reduced when the parent and grand-daughter nuclei 
have different deformations [H, [l3| . 



We have chosen to study the (unphysical) transition between the mirror nuclei ^^Ge and ^^Se in order to have a 
clearer view of the effect of the deformation in the NME's. This transition has the peculiarity that the wave functions 
of the initial and final nuclei are identical (provided Coulomb effects are neglected) and consequently it is easier to 
disentangle the contributions of the Qvf5l3 operator and the nuclear wave functions to the NME. The calculations are 
carried out in the valence space r3g with the effective interaction gcn28:50. The SRC are modeled by a Jastrow factor 
with the Spencer and Miller parametrization [l^ . although it has been shown recently that, once the finite size of the 
nucleon has been taken into account by a dipole form factor, softer options are more realistic [isl [l6| 
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FIG. 1. "^Ge — > ^'^Se NME, M'"", as a function of the difference in deformation induced by the extra quadrupole interaction 
added to '"^Se. 

To increase the deformation of a given nucleus we add to the effective interaction a term X Q ■ Q. Fig. [1] shows 
the results when the final nucleus has been artificially deformed by adding an extra quadrupole-quadrupole term. 
Notice in the first place that for A=0 both nuclei are deformed with /3 ~ 0.2. In spite of that, the NME is a factor 
of two larger than the values obtained for the A=76 and A=82 decays in the same valence space and with the same 
interaction. Hence, even if the two A=66 partners are deformed, the fact that their wave functions are identical 
enhances the decay. Nevertheless, the NME is still far from its expected value in the superfluid limit (NME~8). The 
figure shows that the reduction of the NME as the difference in deformation increases is very pronounced. For the 
values of A between 0.0 and 0.2, the difference in deformation parameter between parent and grand daughter grows 
from zero to about 0.1. In addition, the NME follows closely the overlap between the wave function of one nucleus 
obtained with A=0 and the wave function of the same nucleus obtained with A^^O. This means that, if we write the 
final wave function as: j^*) = aj^'o) + b\'^qq), the 0i^/3/3 operator does not connect and ^gg. This behavior of 
the NME's with respect to the difference of deformation between parent and grand daughter is common to all the 
transitions between mirror nuclei that we have studied (A=50, A=110) and to more realistic cases like the A=82 
decay that we have examined in detail in [l7j . Therefore we can submit that this is a robust result. Similar results 
hold also for the 2z/ decays. 



III. THE NME'S AND THE SENIORITY STRUCTURE OF THE NUCLEAR WAVE FUNCTIONS 

We can also analyze the results of the preceding section in terms of the seniority structure of the wave functions 
of parent and grand daughter nuclei. Indeed when A/3=0 both ^^Ge and ^^Se have identical wave functions. The 
probabilities of the components of different seniority are given in table HI It is seen that changing /3 from 0.22 
(mildly deformed) to 0.30 (strongly deformed) increases drastically the amount of high seniority components in the 
wave function, provoking a seniority mismatch between the decaying and the final nuclei. This leads to very large 
cancelations of the nuclear matrix elements of the decay, as shown also in table HI 

Coming back to the physically relevant decays, we compare in figure [5] the ISM and QRPA NME's. In both 
approaches, the SRC are taken into account in the UCOM framework and gA=1.25 is adopted. We have discussed 
elsewhere that the discrepancies between both approaches show the following trends: when the nuclei that participate 
in the decay have a low level of quadrupole correlations, as in the decays of ^^Zr, ^^^Sn and ^'^^Xe, the calculations 
tend to agree. On the contrary, when the correlations are large, the QRPA in a spherical basis seems not to be able 
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TABLE I. The 


seniority structure of the wave 


functions 


in the A=66 mirror decay 
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s = 4 
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s = 8 


s = 10 


A/3=0 
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43 
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10 
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A^=0.08 
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32 


21 


31 


10 
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M^" 






A/3=0 


-2.02 


3.95 


0.08 


5.16 




A^=0.08 


-0.76 


1.65 


0.02 


2.12 





ISM: full(squares) , s =4(circles) ; QRPA: Tu(bars) , Jy(diamonds) 
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FIG. 2. The neutrinoless double beta decay nuclear matrix elements M' ^^^^ for ISM and QRPA calculations treating the SRC 
with the UCOM approach. Tu, QRPA results from ref. ^ and Jy, QRPA results from refs. 0,01 ■ The ISM results for A=96 
and A=100 are preliminary 



to capture them fully. As the effect of the correlations is to reduce the NME's, the QRPA produces NME's that 
are too large in ^^Ge, ^^Se, ^"^Mo, ^^®Te, and ^^°Te. Indeed, when the ISM calculations are truncated to maximum 
seniority s„i=4, which is the leading order of the ground state correlations in the QRPA (corresponding to the two 
quasi-particle contribution), they follow closely the QRPA results, as can be seen also in figure [5] Notice that only 
when the ISM calculations are converged at this level of truncation the two approaches do produce similar NME's. 

We compare in table |ll] the seniority structure of the wave functions of the ISM and QRPA, in some of the cases 
for which the latter are available (20j . It is seen that the differences are important and share a common trend: in 
the QRPA, the seniority structure of parents and grand daughters is much more similar than in the ISM. According 
to what we have seen in the A=66 case, this is bound to produce larger NME's in the QRPA than in the ISM, as it 
is actually the case. To make this statement quantitative, we have developed the ISM matrix elements in a basis of 
generalized seniority 



Mf,GT,T = ^A^^(a)B„^li3){iyf{(3)\OF,GT,TWi{a)) 
a,/3 



5 





TABLE II 


The seniority structure of tiie 


wave functions 


in the ISM and QRPA 
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where the A's and B's are the amphtudes of the different seniority components of the wave functions of the initial 
and final nuclei. Obviously, when we plug the ISM amplitudes in this formula, we recover the ISM NME's. But, what 
shall we obtain if we put the QRPA amplitudes instead? Indeed, we get approximately the QRPA NME's! (5.73 
for A=76 and 4.15 for A=82). Therefore as we had anticipated, the seniority mismatch of the initial and final wave 
functions, which is severely underestimated in the QRPA calculations, explains most of the discrepancy between the 
two descriptions. In addition, this result strongly suggests that there is some kind of universal behavior in the NME's 
of the neutrinoless double beta decay when they are computed in a basis of generalized seniority. If this is so, the only 
relevant difference between the different theoretical approaches would reside in the seniority structure of the wave 
functions that they produce. 



TABLE III. The GT NME's of the A=48 decay in the generahzed seniority basis 



48^1 


s = 


s = 4 


s = 6 


.s = 8 


"^Ca s = 


3.95 


-3.68 






48Ca s = 4 


0.00 


-0.26 


0.08 


-0.02 



A very spectacular example of the cancellation of the NME by the seniority mismatch is provided by the "^^Ca 
decay. In Table |lT] we have included also the seniority structures of the two nuclei, and we see that they are very 
different. If we now examine the values of the matrix elements {v f (P)\OGTWiio:)) we find the values listed in Table 
mil There are two large matrix elements one diagonal and another off-diagonal of the same size and opposite sign. If 
the two nuclei were dominated by the seniority zero components one should obtain Mgt~4. If "^^Ti were a bit more 
deformed, Mgt will be essentially zero. The value produced by the KB3 interaction is 0.75 that is more than a factor 
five reduction with respect to the seniority zero limit. Earlier work on double beta decays in a basis of generalized 
seniority (limited to s=0 and s=4 components) showing also this kind of cancellations can be found in ref. plj 
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